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Prediction  of  Total  Electron  Content  Using  the 
International  Reference  Ionosphere 


1.  INTRODUCTION 

The  International  Reference  Ionosphere  (IRI)  is  an  empirical  model  of  the 
ionosphere  based  on  experimental  observations.  It  is  the  outcome  of  work  by  the 
URSI  Working  Group  G4,  which  is  identical  to  the  COSPAR  task  group  on  the  Inter¬ 
national  Reference  Ionosphere.  The  goals  and  status  of  the  IRI  have  been  discussed 
by  Hawer  et  al.  1  The  aim  of  the  IRI  is  to  establish  a  compendium  of  height  profiles 
through  the  ionosphere  for  the  four  main  parameters,  namely  plasma  density, 
temperature  of  ions  and  electrons,  and  ion  composition.  These  parameters  are 
generated  from  a  descriptive  model  containing  reliable  data  that  can  be  used  to 
obtain  average  profiles. 

The  IRI  working  group  is  well  aware  that  the  present  model  is  inadequate  in 
some  areas  and  has  encouraged  tests  of  the  model's  validity.  We  describe  here  a 
test  of  the  model's  ability  to  reproduce  observations  of  total  electron  content  (TEC) 
over  a  wide  range  of  conditions.  If  the  tests  are  successful,  they  will  confirm  the 
validity  of  the  electron  density  profile  N(h),  especially  in  the  region  above  the  peak 
of  the  F2  layer  which  contains  approximately  2/3  of  the  TEC. 


(Received  for  publication  15  September  1983) 

1.  Rawer,  K.  (1978)  Goals  and  status  of  the  International  Reference  Ionosphere, 
Res.  Geophys.  Space  Phys.  16:77. 


Rawe  r“  has  described  the  construction  of  the  1HI  N(h)  profile,  which  is  illus¬ 
trated  in  f  igure  I.  The  key  value  is  the  peak  electron  density  NMJ-2,  which  is 
derived  using  the  CCIR  model.  The  bottomside  profiles  are  mainly  based  upon 
lonograms  reduced  to  true  height  profiles.  The  peak  height  HMi  2  ts  obtained  from 
the  parameter  Ali.'iOOOH  using  an  empirical  relationship  due  to  Hilitza.^  The  top¬ 
side  profile  shape  is  deduced  from  a  summary  of  topside  lonosonde  measurements 
(mainly  Al.OL  ETTL'  as  given  by  Bent  and  Llewellyn.0  However,  because  of  limits 

tions  in  the  presentation  of  these  data,  the  I  K1  uses  a  modified  form  of  the  original 
■> 

data.  “ 


Log  N 


2.  Rawer,  K.  (1981)  Introduction  to  1R1  1979,  In  International  Reference  Ionos- 

sphere  -  IRI  79",  J.  V.  Lincoln  and  R.  ().  Conkright  (Eds.)  Report  LAG-82, 
World  Data  Center  A  for  Solar  Terrestrial  Physics,  Boulder,  Colorado. 

3.  CC1R  U974)  Atlas  of  Ionospheric  Characteristics,  CCIR  Report  340-2, 

International  Union  for  Telecommunications,  Geneva. 

4.  Bilitza,  D.  (  197  7)  Modell  zur  Darstellung  der  Hohe  des  F2  -  Maximums  mit 

Hilfe  des  M-3000  -  Wertes  des  CCIR,  Kleinheubacher  Ber.  20:4. 

-  AV. 

5.  Bent,  R.  B.  ,  and  Llewellyn,  S.  K.  11970'  Description  of  the  1965-71  Ionospheric 

Model  l  sed  in  the  Definitive  Systems  (DC)bS',  DBA-Svstems,  Melbourne,  FIT 


A  test  of  the  1R1,  as  far  as  its  ability  to  calculate  TEC  is  concerned,  has  been 

g 

published  by  McNamara  and  Wilkinson.  These  authors  found  that  at  31°  S,  the  IRI 
yielded  discrepancies  in  TEC  sometimes  exceeding  30%,  but  usually  less  than  20%. 
We  present  here  a  much  more  extensive  analysis  based  on  data  from  15  stations, 
covering  a  range  of  magnetic  dip  angles  14  to  77°. 


2.  METHOD  OF  ANALYSIS 

The  IRI  profiles  are  essentially  monthly  median  profiles  and  are  referenc  '  to 
monthly  median  values  of  foF2  and  M(3000)F2.  These  parameters  are  caicui.  t 
using  the  CC1R  coefficients  for  a  given  month  and  linear  interpolation  in  the  aonth 
smoothed  sunspot  number  R.  ^  Consequently  the  most  appropriate  comparison 
calculated  values  of  TEC  is  with  observed  values  which  have  been  linearly  re  1 
to  R.  However,  this  necessitates  a  long  sequence  of  data  covering  a  solar  c 
and  only  a  few  such  data  sets  are  available. 

Since  it  is  the  IRI  N(h)  profile  which  is  of  most  interest  to  us  here,  we  would 
like  to  test  it  under  conditions  for  which  it  is  the  major  source  of  error.  The  other 
main  potential  source  of  error  is  in  the  calculation  of  foF2,  so  it  is  advantageous  to 
use  observed  values  of  foF2  at  the  sub-ionospheric  point  wherever  possible.  Again, 
the  values  of  foF2  should  be  linearly  related  to  R. 

We  have,  therefore,  adopted  a  three-pronged  approach  in  which  we  first  of  all 
aim  for  the  ideal  test  and  then  successively  lift  the  restrictions.  We  consider  three 
types  of  stations /data; 

111  "A”  Stations  are  those  for  w'hich  a  set  of  TEC  data  and  a  set  of  observed 
foF2  data  are  available  for  a  solar  cycle, 

(2)  "B"  stations  are  those  for  which  corresponding  TEC  and  observed  foF2 
data  are  available  on  a  monthly  basis, 

(3)  "C”  stations  are  those  for  which  the  CCIR  method  must  be  used  to 
calculate  foF2. 

Table  1  lists  all  the  stations  considered,  arranged  in  increasing  order  of  dip 
angle.  For  each  station  we  have  considered  where  possible,  data  for  for  March, 

June,  September,  and  December  at  low  and  high  levels  of  solar  activity.  When  par¬ 
ticular  years  must  be  used,  we  have  used  1976  (R  ~  13)  and  1981  (R  ~  140)  if  the 
corresponding  data  were  available.  The  exact  values  of  K  are  not  important. 


6.  McNamara,  L.  F.  ,  and  Wilkinson,  P.  J.  (1983)  Prediction  of  total  electron  content 
using  the  International  Reference  Ionosphere,  J.  Atmos.  Terr.  Phys.  45:169. 

+  For  convenience  in  this  report,  we  shall  use  the  symbol  R  instead  of  the  usual 
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Table  i.  Geographic  i.  orations  vi  the  Observing  Stations,  Sub  -  iop.osplse ri  _■  points 
and  lonosontle  Stations  for  the  P>  Cases  Considered 


Station  |  i.  Coin! 


l.t  Sb 

'  '  T 

i  ■  IiL'  - 

1  Ji’.. 

1  .uni;  . 

r~ 

\  • 

J  ..Tit  • 

ii. 

Manila 

5**. 

121.  ! 

.  . 

.  :  .  H 

12L3 

-L., ml  . 

.  i  . 

121.  ; 

r. 

Ascension 

Island 

-»-o 

345.6 

30.  4 

-7.4 

345.  5 

j 

B. 

Lunping 

25.  2 

122.  2 

35.  2 

23.  0 

121.8 

Chung  Li 

25.  9 

121.2, 

B. 

Palehua 

20.  7 

203.  7 

38.  7 

19.  8 

202.  6 

Maui 

20.  !i 

203.  5 

C. 

Os  an 

37.  1 

127.  0 

5  1.0 

34.  0 

127.  4 

c. 

Ramey 

18.  5 

292.  8 

51.3 

17.  1 

292.  6 

C. 

Athens 

23.  6 

54. 5 

34.  1 

18.  4 

c. 

La  Posta 

32.  5 

242.  9 

58.  0 

30.  0 

245.  0 

c. 

Patrick 

28.  2 

279.  4 

60.  0 

26.  3 

280.  4 

C. 

Shemya 

52.  5 

174.  1 

63.  0 

46.  8 

181.7 

c. 

Boulder 

40.  0 

254.  7 

68.  3 

36.  3 

259.  6 

A. 

Hamilton 

42.  6 

289.  2 

12.  9 

38.  7 

289.  3 

W  allops  1. 

31.9 

284.  5 

c. 

Anchorage 

61.2 

210.  2 

73.  7 

52.  0 

224.  0 

A. 

Goose  Bay 

53.  5 

299.  5 

76.  9 

•4  7 .  5 

297.  8 

St.  Johns 

47.  6 

307.  3 

B. 

N'a  rssa  rssuaq 

61.2 

314.  6 

7  7.5 

53. 1 

307.  8 

Goose  Bay 

53.  3 

299.  5 

The  1  HI  has  been  made  available  as  a  l-'OUTKAN  computer  program  by  World 
Data  Center  A  for  Solar  Terrestrial  Physics.  1  W'e  have  modified  tile  WDC-A  ver- 

g 

sion  to  permit  the  optional  use  of  the  Chiu  model  of  fol'2  and  HMF2  and  of  the 

Cl 

original  Bent  topside  profile  above  HMF2.  "  Listings  of  the  Chiu  model  are  given 
by  Lincoln  and  Conkright.  '  The  modifications  to  include  the  original  Bent  topside 
profile  were  taken  from  the  program  SLANTEc'  with  some  extra  coding  required 
because  of  different  systems  of  units. 

It  would  be  expected  that,  in  general,  the  1K1  would  yield  values  of  TEC  some¬ 
what  less  than  the  observed  values,  since  the  1K1  sets  the  electron  density  to  zero 
above  1000  km.  The  TEC  observations  were  obtained  using  the  Faraday  rotation 
technique,  which  yields  the  TEC  out  to  about  2000  km.  11  Tests  using  the  Bent 
ionospheric  model  indicate  that  the  altitude  range  1000  to  2000  km  contributes 
<  5^0  of  the  TEC  up  to  2000  kin.  This  discrepancy  would  not  be  noticeable  in  the 
present  context. 


(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  43.  ) 
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3.  “A"  STATIONS 


These  are  the  stations  for  which  a  long  sequence  of  TEC  data  covering  a  solar 
cycle  is  available,  in  conjunction  with  a  sequence  of  fol-  2  data  at  the  sub-ionospheric 
point,  also  covering  a  solar  cycle.  There  are  in  fact  only  two  such  stations, 

Hamilton  and  Goose  Bay. 

For  each  of  these  stations,  we  have  correlated  the  observed  monthly  mean 
values  of  TEC  at  each  hour  with  the  smoothed  sunspot  number  R.  The  least  squares 
fit  regression  lines  have  then  been  used  to  determine  the  values  of  TEC  corres¬ 
ponding  to  II  -  0  and  It  ■  100.  A  similar  technique  was  applied  to  foK2  at  the  sub- 
ionospheric  point. 

It  was  found  that  the  values  of  TEC  for  K  -  0  were  unreliable  at  night  when  the 
values  were  very  low.  Consequently,  we  have  been  forced  to  use  actual  observa¬ 
tions  for  low  solar  activity  in  lieu  of  R  -  0  results. 

:f|  llaiiiilloii/Wallnps  Island 

The  Hamilton  l'EC  that  was  used  covered  the  period  1967  to  1980.  The  sub- 
ionospheric  point  lies  very  close  to  Wallops  Island  (see  Table  11  for  which  foE2 
covering  the  period  1975  to  1981  were  used.  The  correlation  coefficients  of  TEC 
and  fol'2  with  R  were  usually  -  0.  95. 

Figures  2  and  1  illustrate  the'  observed  and  calculated  diurnal  variations  of 
TEC  for  Hamilton  at  solar  minimum  (  1976,  R  -  12'  and  maximum  (R  -  100  exactly), 
for  March,  June,  September  and  December.  ^ 

It  can  be  seen  that  there  are  no  major  discrepancies  between  the  observed  values 

of  TEC  and  those  calculated  using  the  IR1.  |Note  that  for  1976,  Hamilton  must  be 

considered  a  "B"  station.  ]  In  general,  the  discrepancies  for  R  =  100  are  least 

when  the  R  -  100  values  of  fol-  2  at  Wallops  Island  are  used.  They  are  greatest  in 

December,  R  100,  between  about  08-12  LT.  The  figures  also  include  the  results 
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for  a  profile  in  which  the  IR1  topside  profile  is  replaced  bv  the  Bent  profile. 


^In  all  diagrams,  the  TEC  is  plotted  as  a  function  of  local  time  in  units  of 
15  -  2 

10  el  m  “,  and  the  following  code  is  used  for  the  data  points; 
x  Observed  monthly  mean  value  (OBS) 

•  Value  given  by  the  original  IRI  (1R1) 

o  Value  given  by  the  IRI  bottomside  and  Bent  topside  models.  This  may  be 
the  value  obtained  using  the  CCIK  value  of  foE2  (IRI-BENT)  or  using  the 
observed  value  of  foE2  at  the  sub-ionospheric  point  (1R1-BENT-OBS) 

+  Value  given  by  the  IRI  using  the  observed  value  of  foE2  (IR1-OBS) 

Where  necessary,  the  data  points  for  the  predicted  values  have  been  slightly  offset 
from  the  corresponding  hour. 
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3.2  Goose  Bay/St.  Johns 

The  Goose  Bay  TEC  data  covered  the  period  1973  to  1980,  while  the  St.  Johns' 
foF2  data  covered  the  period  1970  to  1978.  Figures  4  and  5  show  the  calculated  and 
observed  diurnal  variations  of  TEC  for  Goose  Bay  for  solar  minimum  (R  ~  13)  and 
solar  maximum  (R  =  100)  respectively.  The  greatest  discrepancies  occur  for 
December,  R  =  100,  between  about  08  and  12  LT  (as  with  Hamilton).  The  accuracy 
of  the  CCIR  predictions  of  foF2  for  St.  Johns  is  such  that  the  use  of  "observed" 

R  =  100  values  of  foF2  does  not  lead  to  a  significant  reduction  in  the  discrepancies 
(Figure  5). 

For  the  low  solar  activity  data.  Goose  Bay  may  be  considered  to  be  a  "B"  station 
and  comparisons  made  using  the  monthly  median  values  of  foF2  observed  at  St.  Johns. 
Figure  4  shows  that  when  this  is  done,  the  discrepancies  between  the  observed  and 
predicted  values  tend  to  increase.  This  is  possibly  because  St.  Johns  is  not  quite 
at  the  sub-ionospheric  point  for  Goose  Bay,  but  it  is  also  possible  that  the  better 
results  given  using  the  predicted  values  of  foF2  are  the  result  of  fortuitous  cancella¬ 
tion  of  errors  in  foF2  and  in  the  N(h)  profile. 


4.  “B”  STATIONS 


"B"  stations  are  those  for  which  corresponding  values  of  TEC  and  foF2  at  the 
sub-ionospheric  point  are  available,  but  the  data  do  not  cover  a  solar  cycle.  There 
are  four  such  stations,  Manila,  L. unping  (Taiwain),  Palehua,  and  Narssarssuaq, 
which  have  sub-ionospheric  points  near  the  ionospheric  stations  at  Manila,  Chung  Li, 
and  Goose  Bay  (see  Table  1).  We  consider  Narssarssuaq  first  because  it  causes  the 
least  trouble. 

4.1  Nars*«rMU«q /Goose  Bay 

Narssarssuaq  data  are  available  for  only  1972-  1974,  so  we  consider  data  for 
"low"  activity  (1974,  R  ~  30,  Figure  6)  and  "medium"  activity  (1972,  R  -  ~60, 

Figure  7).  Figure  6  and  7  show  that  the  values  of  TEC  given  by  the  JRI  are  quite 
accurate  in  March,  April,  June  and  September,  especially  when  the  observed 
Goose  Bay  foF2  values  are  introduced.  There  is  a  tendency  for  the  diurnal  varia¬ 
tions  of  the  TEC  based  on  the  observed  values  of  foF2  to  be  somewhat  irregular, 
corresponding  to  irregular  variations  in  foF2  due  to  the  low  sample  sizes  and  the 
complexities  of  the  Goose  Bay  ionograms.  A  discrepancy  between  the  observed 
and  calculated  values  of  TEC  remains  for  December  in  both  1972  and  1974,  between 
08  and  14  L.T,  even  when  the  observed  foF2  values  are  used.  This  indicates  an 
error  in  the  N(h)  profile  given  by  the  IRI. 
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Figure  5.  Predicted  and  Observe  i  Values  of  TEC  at  Goose  Bay  for  Four  Months 
With  High  Solar  Activity  (R  =  100) 
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Figure  7.  Predicted  and  Observed  Values  of  1  PC  at  Narssarssuaq  for  Four 
Months  With  Medium  Solar  Activity  (1972) 
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4.2  Manila/Manila 


It  is  in  the  calculation  of  the  Manila  TEC  that  we  encounter  the  first  major  dis¬ 
crepancies  between  the  calculated  and  observed  values  of  TEC.  Manila  TEC  data 
were  available  for  1980-1982,  apart  from  missing  months.  Figure  8  shows  the 
observed  and  predicted  values  of  TEC  for  April  1982  (H  =  124),  June  1981  (H  =  155), 
October  1981  (R  =  142),  and  December  1981  (R  =  138).  (March  and  September  data 
were  not  available. ) 

In  general,  the  IRI  predicts  daytime  values  of  TEC  which  are  only  about  one- 
half  of  the  observed  values,  in  all  months.  The  errors  at  night  are  somewhat 
smaller.  The  source  of  the  error  does  not  lie  in  the  prediction  of  fol'2,  since  as 

Figure  9  shows,  the  error  in  the  predicted  foF2  is  typically  only  about  10%.  Use 

9 

of  the  original  Bent  profiles  improves  the  situation  but  still  leads  to  prediction 

15-2 

errors  reaching  300  x  10  el  m  and  incorrect  diurnal  variations. 

4.3  Limping/ Chung  Li 

The  Lunping  TEC  data  were  available  for  1980-March  1983,  with  some  missing 
months.  We  have  compared  the  predicted  and  observed  values  of  TEC  for  March 
1980  (K  -  80),  June  1981  (R  =  142),  September  1982  (R  =  109),  and  December  1981 
1R  =  138). 

Figure  10  shows  that,  as  with  Manila,  the  1R1  greatly  underestimates  the  TEC 
in  some  months,  that  the  prediction  of  fol'2  is  not  the  problem,  and  that  the  use  of 
the  Bent  topside  profile  leads  to  more  accurate  results.  However,  the  discrepancy 
between  the  predicted  and  observed  TEC  values  still  exceeds  300  >  lO1'’  el  m  ^ 
during  the  day  in  September  and  December. 

4.4  Paleliua/Maiii 

The  l’alehua  TEC  data  were  available  for  1980-1982  and  we  have  compared  the 

predicted  and  observed  TEC  values  for  March,  June,  September,  and  December 

(K  -  140).  Figure  II  shows  that  l’alehua  follows  i. unping  and  Manila  in  the  way  the 

IRI  greatly  underestimates  the  TEC  during  the  day  during  some  months.  The  dis- 

15  -2 

crepancv  is  least  for  June,  when  it  is  of  the  order  of  80  >  10  el  m  or  about  15%. 
The  daytime  discrepancies  remain  essentially  unchanged  when  Maui  observations 
of  foE2  are  used.  However,  the  discrepancies  at  night  tend  to  decrease. 

Use  of  the  original  Bent  profiles  above  the  peak  of  the  layer  yields  values  in 
better  agreement  with  the  observations,  but  large  discrepancies  still  remain  in 
December.  The  calculated  diurnal  variation  agrees  well  with  the  observations  in 
September,  but  not  in  the  other  months. 
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Figuro  11.  Predicted  am)  Observed  Values  of  TKC  at  Palehua  (Hawaii)  for  Four 
Months  With  High  Solar  Activity  (1981) 


5.  "C”  STATIONS 


These  are  the  stations  for  which  the  calculation  of  foE2  relies  on  the  CC1R 
method,  with  no  real  means  of  checking  the  calculated  values.  However,  the  re¬ 
sults  for  the  "A"  and  "B"  stations  suggests  that  the  major  source  of  any  discrepan¬ 
cies  will  lie  in  the  topside  N(h)  profile.  We  consider  the  "C"  stations  in  order  of 
increasing  dip  latitude. 

5. 1  Ascension  Island 

For  Ascension  Island  we  have  considered  the  months  March  1981  <K  -  143), 

June  1980  (R  =  153),  September  1980  (R  =  150),  and  December  1980  (R  =  143).  In 
general,  the  daytime  TEC  values  predicted  by  the  I  HI  are  too  low  by  at  least 
200  x  lO1,0  el  m  ,  the  discrepancy  reaching  800  v  101"1  el  m  “,  or  '  50"'[>,  in 
March  1981  (see  Figure  12).  The  original  Bent  topside  profile  yields  somewhat 
more  accurate  results,  but  for  June  1980,  retains  the  peculiar  diurnal  variation 
given  by  the  I  HI. 

The  post-sunset  peak  in  the  observed  values  of  TEC  is  missing  from  both  sets 
of  predicted  values,  although  there  is  a  plateau  from  1800-2400  l.T  in  December  1980. 

5.2  Osan.  hurra 

TEC  data  for  Osan  were  available  for  1980-1982,  with  some  missing  months, 
noticeably  the  equinoxes.  Wo  have  considered  the  data  for  April,  June,  October 
and  December  1“K1  <R  -  140)  and  the  results  are  shown  in  f  igure  13. 

The  daytime  values  of  TEC  calculated  using  the  I  111  are  generally  quite  good, 

i  y  -  > 

although  the  errors  reach  ~  140  s  10  '  el  m  “  for  October  1981.  The  nighttime 

values  given  bv  the  1R1,  on  the  other  hand,  tend  to  be  too  high  bv  a  factor  of  ~  2, 

1 5  -  2 

with  a  typical  absolute  error  of  ~100  ^  10  el  m 

5.3  Rainer 

Ramey  TEC  observations  wore  available  for  1980-1982  and  we  have  considered 

the  data  for  March,  June,  September,  and  December  1981  t  1  {  -  1401,  J'he  predicted 

values,  shown  in  f  igure  14,  tend  to  underestimate  the  observed  TfX  by  up  to 

'  200  ■  10  ’’  el  ni  “  i  ~  2.V'i>)  during  the  day,  except  in  June.  There  is  no  consistent 

discrepancy  at  night,  with  only  a  small  discrepancy  in  June  and  December,  but  an 

15  -  2 

underestimate  bv  -  80  •  10  el  m  «~30"»)  during  March  and  September. 
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Athens  TEC  data  were  available  for  1980-1982,  with  some  missing  months. 

We  have  considered  the  data  for  March  1982  (R  =  129),  June  1981  (R  =  142), 
September  1982  (R  =  117),  and  December  1982  (R  =  138).  The  results  are  shown  in 
Figure  15. 

The  1R1  TEC  values  during  the  day  tend  to  be  too  low  by  up  to  ~  200  X  lO1^  el  m” 

in  March  and  December,  while  the  values  during  the  night  tend  to  be  too  high  by  up 
15-2 

to  ~  100  x  10  el  m  in  September  and  December. 

5.5  La  Post  a 

La  Posta  TEC  data  were  available  from  1980  to  1983  and  we  have  considered 
the  four  months  March,  June,  September,  and  December  1981  (R  ~  140)  (see  Fig¬ 
ure  16).  During  the  day,  the  discrepancies  between  observed  and  predicted  values 
15-2 

reach  -  150  X  10  el  m  in  March  and  December,  with  the  1R1  values  being  too 
low.  At  night,  the  predicted  TEC  values  tend  to  be  too  high  by  a  factor  '  2.  Fig¬ 
ure  16  also  shows  the  TEC  values  calculated  using  the  original  Bent  topside  profiles. 
In  this  case,  the  results  are  very  similar  to  the  1R1  results. 


5.6  Patrick 

TEC  data  for  Patrick  were  available  for  1980-1983  and  we  have  considered 

March,  June,  September,  and  December  1981  (R  -140)  (see  Kigure  17).  In  general, 
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the  IRI  underestimates  the  TEC  during  the  day  by  up  to  '200  >10  el  m  except 
in  June  when  the  discrepancies  are  very  low.  The  discrepancies  at  night  do  not 
yield  a  consistent  picture. 


5.7  Shetnya 

For  Sheyma,  TEC  data  were  available  for  September  1977  to  April  197°,  and 
for  1980-1982.  We  have  considered  the  months  March  1979  (R  =  137),  June  1978 
(R  =  89),  September  1978  <K=  108),  and  December  1978  <H  =  118)  (see  Figure  18), 

In  general,  the  discrepancies  are  reasonably  small,  except  in  June  1978,  when  the 
calculated  and  observed  diurnal  variations  of  TEC  do  not  follow  each  other  very 
closely. 

5.8  Boulder 

Boulder  TEC  data  were  available  for  1980-1982  and  we  have  considered  the 
months  March,  June,  September,  and  December  1981  (R~14Q).  In  general,  the 
agreement  between  the  observed  and  predicted  values  of  TEC  is  quite  good  during 
the  day,  but  the  predicted  nighttime  values  are  a  factor  of  ~2  too  high  (see  Figure  19). 
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Figure  IP.  Predicted  and  Observed  \  alues  of  TEC  at  Boulder  for  Four  Months 
With  High  Solar  Activity  (1980 
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5.9  Anchorage 

Anchorage  TEC  data  were  available  for  1980-1983,  with  some  missing  months, 
and  we  have  considered  March  1981  (R  =  143),  June  1982  (R  =  117),  September  1981 
(R  =  143),  and  December  1980  (R  =  143)  (see  Figure  20).  The  discrepancies  between 
calculated  and  observed  values  of  TEC  during  the  day  show  no  consistent  trend,  the 
IRI  values  being  too  low  in  March,  too  high  in  December  and  fairly  minor  in  June 
and  September.  The  discrepancies  at  night  are  also  fairly  minor.  In  this  case, 
the  IRI  values  match  the  observations  better  than  do  the  values  calculated  using  the 
original  Bent  topside  profile. 


6.  DISCUSSION  OF  RESULTS 

Since  we  have  found  that  errors  in  the  calculated  values  of  foF2  are  not  the 
major  cause  of  the  discrepancies  between  calculated  and  observed  values  of  TEC, 
and  that  the  discrepancies  are  significantly  greater  at  low  latitudes,  we  shall  discuss 
the  results  under  three  latitude  headings-mid,  high,  and  low. 

6.1  Mid-Latitude  Stations 

We  consider  here  Osan,  Ramey,  Athens,  La  Posta,  Patrick,  Shemya,  Boulder, 
and  Hamilton.  The  low  latitude  limit  was  set  by  the  fact  that  Palehua  TEC  is 
affected  by  the  equatorial  anomaly,  while  the  upper  latitude  limit  was  set  to  L  =  4. 

The  Boulder,  Hamilton,  and  Osan  results  show  no  major  discrepancies  between 
the  observed  and  calculated  values  of  TEC  during  the  day,  but  the  IRI  values  tend 
to  be  too  high  at  night  by  a  factor  '  2.  The  IRI  results  for  Athens,  La  Posta, 

Patrick,  and  Ramey  are  too  low  during  the  day  in  March,  September,  and  December, 
but  have  only  minor  errors  in  June  (summerl.  There  is  no  consistent  discrepancy 
at  night.  The  discrepancies  at  Shemya  are  similar,  except  for  June,  when  the  cal¬ 
culated  values  of  TEC  vary  substantially  both  in  absolute  value  and  diurnal  variation 
from  the  observed  values.  The  cause  of  these  discrepancies  is  not  known,  but  the 
double-peaked  diurnal  variation  of  the  TEC  (which  is  also  found  for  June  19811 
indicates  that  in  June  the  ionosphere  at  Shemya  differs  from  the  normal  mid-latitude 
ionosphere.  The  discrepancies  at  the  other  stations  in  the  other  seasons  may  in  fact 
be  partly  due  to  errors  in  the  predicted  values  of  foF2,  but  the  consistency  of  the 
errors  over  a  wide  range  of  longitudes  also  suggests  that  the  profile  shape  may  be 
at  fault  as  well. 
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6.2  High-ljlitude  Stations 


The  results  for  Anchorage,  Goose  Bay,  and  Narssarssuaq  are  similar  to  those 
for  the  mid-latitude  stations,  with  no  major  discrepancies,  especially  when  ob¬ 
served  values  of  foF2  are  available.  The  largest  discrepancy  occurs  for  about 
08-14  LT  in  December  when  the  I  HI  values  are  too  high. 

6.3  l>ow- Latitude  Stations 

The  four  low  latitude  stations,  Manila,  Ascension  Island,  1. unping,  and  Falehua 
all  yield  considerable  discrepancies  between  the  observed  and  calculated  values  of 
TEC.  During  the  day,  the  calculated  values  can  be  a  factor  of  2  too  low  (absolute 
error  <  500  >  10  '  el  m  >,  with  an  incorrect  diurnal  variation.  The  1HJ  also 
does  not  reproduce  the  observed  post-sunset  increase  of  TEC  at  Ascension  Island. 

As  indicated  by  the  results  for  Manila,  1. unping,  and  l’alehua,  the  discrepancies 
are  not  due  to  errors  in  the  predicted  values  of  fol  2  and  must  therefore  be  due  to 
the  N(h)  profile  shape.  I.arger  daytime  values  of  TEC,  which  agree  better  with  the 
observed  values,  are  obtained  by  replacing  the  I  HI  topside  profiles  bv  the  original 
Bent  profile.  While  this  procedure  was  not  found  necessary  at  mid-  and  high- latitude 
stations,  it  is  at  least  a  step  in  the  right  direction  for  low  latitude  stations.  In  some 
months,  the  modified  I  HI  profiles  yield  good  agreement  with  the  observations.  Even 
so,  the  calculated  diurnal  variations  of  TEC  often  do  not  match  the  observed  day¬ 
time  variations,  nor  the  post-sunset  increases  at  Ascension  Island. 

The  IRI-BENT  (that  is,  with  the  original  Bent  topside  profile)  results  for  Manila 
indicate  a  "crisis  of  confidence"  which  occurs  at  around  08  LT— the  calculated  TEC 
flattens  out  rather  than  continuing  the  post-sunrise  increase.  Consideration  of  the 
diurnal  variation  of  foE2  at  Manila,  with  foF2  going  through  10  MHz  at  around  07  LT, 
led  us  to  consider  a  possible  modification  of  the  Bent  topside  profile,  in  particular 
the  semi-thickness  of  the  topside  parabolic  YT.  The  coding  for  FKOFL2  (Llewellyn 
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and  Bent;  especially  page  181)  shows  that  if  foF2  exceeds  10.5  MHz,  YT  is  given 
by  the  formula: 

YT  =  YM  X  {0.  133  X(foF2  -  10.  5)  +  l}  , 

where  YM  is  the  half  thickness  of  the  bottom«ide  parabola.  Since  the  TEC  increases 
when  YT  is  increased,  the  value  of  TEC  given  by  the  IRI-BENT  model  can  be  in¬ 
creased  when  foF2  >  10.5  MHz  simply  by  increasing  the  factor  0.  133.  Increasing 
it  to  0.  5  yields  TEC  values  more  in  agreement  with  the  observed  daytime  values. 

The  predicted  daytime  TEC  values  for  the  IRI-BENT  model  can  also  be  in¬ 
creased  (by  a  lesser  amount)  by  relaxing  the  upper  limit  on  the  10.  7  cm  flux  set  by 
Llewellyn  and  Bent.®  Subroutine  REFRAC  (Report  page  172)  limits  the  flux  to  130. 
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Removing  this  limit  allows  the  exponent  of  the  lower  topside  exponential  to  achieve 
lower  values,  yielding  a  more  slowly  decaying  profile.  (See  Figure  4  of  Llewellyn 
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and  Hent.  )  The  exponents  must,  however,  be  constrained  to  stay  positive,  which 
w..s  the  original  reason  for  the  limit  on  the  flux  value. 

Hotter  agreement  with  the  observed  daytime  TEC  values  at  low  latitude  stations 
may  thus  he  obtained  by  judicious  manipulation  of  the  parameters  of  the  Bent  top¬ 
side  profile.  We  have  not  pursued  this  point.  All  results  given  in  this  report  were 
obtained  using  the  unmodified  Hent  profiles. 

The  rill  at  Ascension  Island  has  been  modeled  theoretically  bv  Anderson  and 
1 2 

Klobuchar,  who  were  able  to  reproduce  the  observed  TEC  variations  including 
the  late  afternoon  decrease  and  post-sunset  enhancements.  In  an  (as  vet)  incomplete 
study  of  the  effects  of  A  _U  drifts  at  equatorial  stations,  Anderson  and  McNamara 
calculated  the  Nth)  profiles  and  TKC  for  Palmyra  (5.  9  N,  197.  8  E,  13°  dip)  for 
January  1949  (R  =  137).  The  calculated  TKC  values  are  shown  in  Figure  31,  which 
also  shows  the  observed  TKC  for  Manila  for  January  1982  (R  =  137),  the  I  HI  TEC 
values  and  the  daytime  IR1-BENT  TKC  values.  The  theoretical  values  (A)  are 
clearly  the  best  match  to  the  observed  data,  although  some  discrepancies  still 
remain  at  about  18-22  I  T.  These  may  in  fact  be  due  to  different  17  *  H  drifts 
existing  at  Manila  and  Palmyra,  and  are  not  our  present  concern.  We  are  con¬ 
cerned  rather  with  the  considerably  different  N(h)  profiles  given  by  the  1RI-BKNT 
program  and  by  the  theoretical  model.  Figure  22  shows  that  the  two  profiles  for 
12  I  l',  January,  R  -  137,  bear  little  resemblance  to  each  other,  the  theoretical 
profile  being  much  broader  than  the  1RI-BENT  profile  (and  consequently  yielding  a 
greater  TKC).  This  discrepancy  points  up  the  need  for  further  studies  of  the  pro¬ 
file  shape  at  low  latitudes.  The  fact  that  at  least  half  of  the  "missing''  TKC  would 
appear  to  come  from  below  the  peak  of  the  layer  suggests  caution  in  trying  to  match 
the  observe'  dues  of  TKC  simply  by  adjusting  parameters  of  the  topside  profile. 

6.4  Nighttime  Diurrepaneie* 

We  have  tended  so  far  to  discuss  mainly  the  discrepancies  which  occur  during 
the  day,  since  these  have  been  found  to  have  the  largest  absolute  values.  However, 
the  discrepancies  might  tend  to  reach  higher  relative  values,  which  may  be  of  more 
concern  in  some  applications  than  the  absolute  discrepancies.  There  seems  to  be 
no  general  trend  in  the  discrepancies,  overestimates  and  underestimates  being 
about  equally  likely.  Some  of  the  "A"  and  "B"  stations,  (for  example,  Hamilton  and 

12.  Anderson,  D.N.,  and  Klobuchar,  J.A.  ( 1983)  Modeling  the  total  electron 

electron  content  observations  above  Ascension  Island,  J.  Geophys.  Res.  , 
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(.loose  Bay)  shows  very  small  discrepancies  when  observed  values  of  fol'2  are  used. 
For  Manila  and  Narssarssuaq,  on  the  other  hand,  the  IKI  systematically  under¬ 
estimates  the  TEC.  even  when  the  observed  foF'2  values  are  used. 


TEC 


Figure  21.  Calculated  and  Observed  Values  of  TEC  lor 
January  at  a  Station  at  13°  Dip,  for  H  =  137.  The  ob¬ 
served  values  are  for  Manila,  while  calculated  values 
are  given  for  the  IKI,  the  modified  I  HI- BENT  and  the 
Anderson  theoretical  models.  The  Anderson  values 
were  calculated  for  Palmyra 
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7.  CONCLUSIONS 


We  have  found  that,  in  general,  the  I  HI  yields  reasonable,  and  often  quite 
accurate,  values  of  TEC  at  mid-  and  high-latitude  stations,  but  that  it  tends  to 
underestimate  severely  the  daytime  values  of  TEC  at  low  latitude  stations.  Better 
agreement  with  the  observed  daytime  values  of  TEC  at  low  latitudes  can  be  obtained 
by  replacing  the  IRI  topside  profile  by  the  original  Bent  profile,  but  theoretical 
model  studies  indicate  that  a  substantial  revision  of  both  tne  1  HI  and  I  HI- BENT 
profiles  is  required  at  low  latitudes.  In  those  areas  for  which  agreement  has  been 
found  between  observed  and  predicted  values  of  TEC,  the  I  HI  electron  density  pro¬ 
files  are  presumably  good  representations  of  the  actual  profiles.  From  an  opera¬ 
tional  point  of  view,  this  means  that  the  I  HI  may  be  used  to  calculate  reliable  values 
of  TEC  in  those  same  areas. 
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